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1 12b.  DISTRBUTION 


Representation  of  natural  terrain  growth  as  a  complex  dielectric  slab  with  some  arbitrary  thickness 
equivalent  to  the  mean  height  of  the  vegetation  is  a  common  modeling  simplification  technique. 
The  vegetation  medium  may  be  viewed  as  two  or  more  different  types  of  dielectric  scatterers 
arranged  in  a  random  geometric  configuration  in  a  surrounding  “host”  medium.  Based  on  this 
assumption,  the  foliage  is  treated  as  a  continuous  medium,  with  the  propagation  of  the  waves 
governed  by  an  effective  dielectric  constant,  (e*).  In  this  paper  several  volumetric  multiple¬ 
scattering  models  were  used  to  determine  the  effective  dielectric  constants  representative  of  foliage 
media.  A  brief  discussion  of  the  Average  T-Matrix  Approximation  (ATA)  and  the  Coherent 
Potential  Approximation  (CPA)  is  presented.  The  ATA  represents  the  field  at  an  individual 
scattering  center  in  terms  of  an  equivalent  field  produced  by  the  other  scatterers  but  includes  no 
provision  for  the  interaction  between  scatterers.  The  CPA  neglects  the  difference  between  the 
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external  excitation  field  and  the  average  field  and  includes  interaction  effects  between  scatterers. 
The  effective  parameters  obtained  using  these  approximations  are  contrasted  with  results  obtained 
using  several  other  mixing  formulas.  In  this  paper  the  multiple  scattering  approximation  techniques 
are  used  to  quantitatively  estimate  e *  for  vegetation  media  based  on  the  foliage  density  and 
dielectric  properties  of  the  scattering  elements  comprising  the  foliage.  Additionally,  a  formalism  is 
presented  that  allows  incorporation  of  the  propagation  losses  experienced  by  an  electromagnetic 
wave  propagating  through  vegetation  media  into  electromagnetic  scattering  models,  and  yet  is  not 
computationally  intensive.  The  theoretical  approximations  of  (*  are  compared  with  data  at  four 
levels.  First,  the  calculated  values  of  e*  were  compared  with  values  reported  in  the  literature. 
Second,  skin  depth  penetration  distances  were  calculated  and  compared  to  experimental  values. 
Thirdly,  the  effective  dielectric  constants  for  a  foliage  environment  were  used  to  calculate 
attenuation  coefficients  of  coherent  waves  propagating  through  dense  vegetation.  The  calculated 
attenuation  constants  were  compared  with  experimental  measurements  reported  by  various  authors 
at  frequencies  between  50MHz  and  3.2GHz,  resulting  in  good  agreement.  Finally,  the  values  of  e* 
were  incorporated  into  a  bistatic  scattering  model  which  was  used  to  calculate  an  effective 
normalized  scattering  cross  section  <r*  for  a  grass  and  forest  covered  terrain.  This  was  compared  with 
L-band  data  resulting  in  excellent  agreement  between  theory  and  experimental  data. 
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1.0  INTRODUCTION 

The  problem  of  determining  the  attenuation  loss  experienced  by  an 
electromagnetic  (EM)  wave  propagating  through  a  foliage  environment  is 
important  and  quite  complex.  One  common  method  used  to  model  this 
phenomenon  is  the  continuous  medium  approach  where  the  attenuation  properties 
in  a  foliage  medium  are  governed  by  an  effective  dielectric  constant  (c*).  A  simple 
model,  applying  dielectric  mixing  formula,  has  been  developed  which  can 
quantitatively  estimate  the  effective  dielectric  constants  of  a  foliage  layer  based  on 
several  descriptive  parameters  of  the  foliage  media.  The  theoretical 
approximations  of  e*  obtained  from  the  model  are  compared  with  data  at  four 
levels.  First,  direct  comparison  of  calculated  effective  parameters  with  values  for 
similar  conditions  reported  in  the  literature.  Second,  values  of  skin  depth 
penetration  calculated  using  the  effective  parameters  were  compared  with 
experimentally  measured  values.  Third,  estimations  of  the  attenuation  loss 
(dB/m)  of  coherent  waves  propagating  through  a  dense  forest  were  compared  with 
experimental  data  reported  by  various  authors  at  frequencies  between  50MHz  and 
3.2GHz.  Finally,  values  of  e*  were  incorporated  into  a  bistatic  scattering  model 
which  was  used  to  calculate  an  effective  normalized  scattering  cross  section  (<r°)  for 
a  sod  field,  grass  and  forest  covered  terrain.  This  was  compared  with  L-band  data 
resulting  in  excellent  agreement  between  theory  and  experimental  data. 

The  model  assumes  that  a  vegetation  medium  may  be  viewed  as  two  or 
more  different  types  of  dielectric  scatterers  arranged  in  a  random  geometric 
configuration  in  a  surrounding  “host”  medium.  The  use  of  an  effective  dielectric 
constant  e*  permits  a  heterogeneous  mixture  of  random  scattering  elements  to  be 
treated  as  a  homogeneous  medium  and  provides  a  simple  mechanism  to  account  for 
multiple  scattering  events  expected  to  occur  when  an  electromagnetic  wave 
propagates  through  that  medium.  This  paper  applies  several  multiple  scattering 
approximation  techniques  to  estimate  values  representative  of  e*  of  a  vegetation 
medium  based  on  readily  measured  quantities  such  as  the  foliage  density  and 
dielectric  properties  of  the  scattering  elements  comprising  the  foliage,  i.e.  wood, 
leaves;  and  presents  a  formalism  that  allows  the  reader  to  incorporate  the 
propagation  loss  of  an  electromagnetic  wave  into  complex  electromagnetic 
scattering  models,  and  yet  is  not  computationally  intensive. 
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A  brief  discussion  of  the  general  theory  of  the  multiple  scattering  models 
that  are  applied  in  determining  the  values  of  c*  is  presented  in  Appendix  A. 
Principal  formula  discussed  include  the  Average  T-Matrix  Approximation  (AT A) 
and  Coherent  Potential  Approximation  (CPA).  The  ATA  represents  the  field  at 
an  individual  scattering  center  in  terms  of  an  equivalent  field  produced  by  the 
other  scatterers  but  includes  no  provision  for  the  interaction  between  scatterers. 
The  CPA  neglects  the  difference  between  the  external  excitation  field  and  the 
average  field  and  includes  interaction  effects  between  scatterers.  In  addition, 
several  other  mixing  formula  are  presented  in  order  to  compare  the  effects  that 
differing  modeling  assumptions  and  variation  of  the  assumed  shapes  of  the 
embedded  scatterers  have  on  the  determination  of  the  effective  dielectric  constants. 

Additionally,  Appendix  B  includes  a  discussion  of  the  methods  used  to 
obtain  the  dielectric  constants  of  the  components  that  comprise  vegetation  media 
(i.e.  leaves  and  wood).  These  values  are  known1,2  to  be  variant  with  respect  to  the 
moisture  content  and  salinity  of  the  wood,  ambient  air  temperature,  type  of  wood, 
and  orientation  of  the  wood  grain  with  respect  to  the  electric  field  polarization. 
Although  variation  of  any  of  these  conditions  would  cause  alteration  of  the 
calculated  effective  dielectric  constants,  detailed  analysis  of  these  effects  are  not 
included  in  this  paper  which  is  intended  to  demonstrate  the  viability  of  using  the 
multiple-scattering  models  to  determine  estimates  of  the  effective  dielectric 
constants  of  vegetation  media,  not  on  determining  the  values  of  £*  for  various 
vegetation  environments. 

2.0  PROPAGATION  OF  EM  WAVES  THROUGH  FOLIAGE 

A  limited  number  of  studies  have  been  done  that  examine  the  multiple 
scattering  and  absorption  events  expected  to  occur  when  electromagnetic  waves 
propagate  in  a  foliage  medium.  Some  treat  each  component  of  the  vegetation  as  an 
individual  scatterer  and  define  a  shape  and  dielectric  constant  associated  with  it1,3,4 
while  others  treat  the  foliage  as  a  heterogeneous  ensemble  of  scatterers  of  arbitrary 
shape  and  orientation5,6,7,8.  In  many  instances,  however,  application  of  these 
processes  in  a  system  simulation  program  is  impractical  due  to  their  complexity  or 
computational  requirements.  Because  of  this,  various  assumptions  are  usually 
introduced.  One  common  simplification  involves  the  representation  of  natural 
terrain  growth  as  a  complex  dielectric  slab  with  some  arbitrary  thickness 
equivalent  to  the  mean  height  of  the  vegetation.  In  this  model  the  vegetation 
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components  are  considered  individual  scatterers  embedded  in  a  host  medium. 
Based  on  this  assumption,  the  foliage  is  treated  as  a  continuous  medium,  with  the 
propagation  of  the  waves  governed  by  an  effective  dielectric  constant,  c*.  When 
considering  propagation  through  and  scattering  from  vegetation  alone,  use  can  be 
made  of  the  experimental  fact  that  the  volume  fraction  occupied  by  the  foliage  is 
usually  less  than  5%,  i.e.  the  foliage  is  sparse9.  Another  assumption  made  in  the 
modeling  of  the  propagation  through  and  scattering  from  foliage  is  that  the 
scattering  from  the  interconnects  between  leaves  and  branches  and  leaves  and 
stalks  can  be  neglected.  As  pointed  out  by  Brown1,  this  is  a  good  assumption  since 
the  component  elements  of  foliage  are  much  better  absorbers  of  EM  energy  than 
scatterers  of  EM  energy  in  the  frequency  range  of  interest  here  (50MHz  to  3GHz). 
Because  the  foliage  is  sparse  and  absorption  rather  than  scattering  is  the  primary 
interaction  process,  the  coherent  field  is  dominant  over  the  incoherently  scattered 
waves  within  the  medium. 

The  total  field  inside  a  vegetation  medium  can  be  mathematically  divided 
into  the  sum  of  a  mean  or  coherent  part,  and  a  zero  mean  fluctuating  part10, 

E  =  <E>  +  SE/juc  ,  (1) 

where  <  -  >  denotes  an  ensemble  average  and  <5E  is  a  zero  mean  field  quantity.  In 
order  to  apply  effective  parameters  to  describe  the  absorption  and  scattering  of 
electromagnetic  waves  from  a  vegetation  environment,  certain  conditions  must  be 
met.  First,  the  average  field  <E>  must  be  much  greater  than  the  fluctuation 
field  SEfiuc  having  zero  mean, 

<E>  >>  SEfiuc  .  (2) 

Secondly,  the  difference  between  the  wave  numbers  of  the  wave  propagating  in  the 
vegetation  medium  and  in  free  space  should  be  less  than  the  reciprocal  of  the 
maximum  dimension  of  the  embedded  scatterers11, 

(k  -  k0)  -  L  <  1.0  .  (3) 

For  spherical  embedded  scatterers  L  will  be  the  on  the  order  of  the  diameter  of  the 
spheres.  The  limitations  of  the  frequency  and  scatterer  dimension  on  the 
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determination  of  the  effective  parameters  of  foliage  media  in  this  study  will  be 
bounded  by, 

(^-l)(^)-L  <  0.5  ,  (4) 

where  t*  is  the  effective  complex  dielectric  constant  of  the  foliage  medium,  L  is  the 
largest  dimension  of  the  embedded  scatterers,  and  A  is  the  wavelength  of  the 
incident  radiation.  The  condition  given  by  Eq.  (4)  essentially  requires  that  the 
scattering  properties  of  a  body  in  free  space  and  in  the  random  medium  not  be 
substantially  different,  a  reasonable  assumption  for  foliage  media,  which  has 
effective  dielectric  constants  close  to  that  of  free  space12. 

If  the  above  conditions  are  met,  e*  can  be  used  to  estimate  the  amplitude  of 
an  electric  field  that  has  propagated  some  distance  through  a  heterogeneous 
random  medium.  The  amplitude  of  the  resulting  electric  field  is  defined  by  the 
relationship, 

E  =  Eq  -exp(jk0>R*-z)  ,  (5) 

—¥  — ♦ 

where  k0=27r/A,  Eq  is  the  amplitude  of  the  incident  electric  field,  E  is  the 

amplitude  of  the  electric  field  after  propagating  through  the  media,  e*  is  the 
effective  complex  dielectric  constant  of  the  media,  and  z  is  the  distance  the  wave 
travels  through  the  media.  The  slab  model  is  most  useful  when  applied  in  a  system 
simulation  program  where  the  foliage  media  represents  an  attenuation  factor  and  is 
not  the  primary  focus  of  the  model.  For  example,  when  modeling  the  scattering 
from  a  vegetation  covered  terrain  the  vegetation  layer  would  be  represented  by  the 
homogeneous  dielectric  slab  that  covers  the  ground,  which  itself  could  be 
represented  by  a  rough  surface  model.  An  example  of  an  application  similair  to  this 
is  discussed  in  the  section  entitled  "Application  of  the  Dielectric  Slab  Model”  later 
in  this  report.  The  dielectric  slab  model  hais  been  shown  to  provide  a  good 
estimate  of  the  scattering  and  absorption  of  electromagnetic  radiation  propagating 
through  foliage  media  when  used  in  such  a  "background  configuration”  up  to  L- 
Band13,14. 

3.0  THEORY 

Various  analytical  techniques  have  been  developed  in  an  attempt  to 
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determine  effective  parameter  values  to  describe  heterogeneous  media  composed  of 
a  host  medium  a r  i  various  embedded  scatterers3’5,15*16.  Few  of  these  methods  have 
been  applied  Lo  the  determination  of  the  effective  parameters  of  foliage  media.  One 
method  that  appears  to  be  well  suited  to  this  application  is  the  recent  work  of 
Blankenship5,  Lang17  and  Kohler  and  Papanicolaou18,  that  extend  the  earlier  work 
of  M.  Lax19  and  is  used  as  the  basis  for  this  analysis. 

Assume  we  are  given  a  region,  or  scattering  medium,  0  comprised  of  free- 
space  with  relative  dielectric  constant  (€0=1)  that  contains  m  different  embedded 
scatterers  with  relative  complex  dielectric  constants  tj  cOT.  Defining  VtJ,  i  = 
l,...m,  j  =  l,2,...Nj,  as  the  subset  of  O  occupied  by  embedded  scatterers  j  of  class  i, 
the  dielectric  properties  of  the  composite  medium  occupying  0  axe  defined  by 
Blankenship2, 


where, 


Ni 

f(x)  =  f0  +  E(Vfo)  £x„y(x) 

1=1  J=1  IJ 


X  €  V.J 
X  *  Vy 


(6) 

(7) 


Assume  that  6,  is  defined  to  be  a  dimensionless  parameter  describing  a  particular 
aspect  of  the  elementary  scatterers  in  class  i,  such  as  the  radius  of  spherical 
scatterers;  and  p  is  the  total  volume  fraction  of  O  occupied  by  scatterers  of  class  i 
=  l,...,m.  Given  a  field  incident  on  the  region  0,  the  objective  is  to  characterize 
the  scattering  properties  of  the  composite  medium  for  the  limit  N,-»  oo,  6  — ♦  0,  and 
where  p.  is  a  constant.  Using  the  transition  operator  convention  to  derive  various 
approximations  it  is  possible  to  solve  for  (*  of  a  foliage  medium  which  can  be 
defined  as  a  host  medium  containing  multiple  embedded  scattering  types.  Two 
such  approximations  are  examined  in  this  paper.  The  Average  T-matrix 
Approximation  (ATA)  which  represents  the  field  at  an  individual  scattering  center 
in  terms  of  an  equivalent  field  produced  by  the  other  scatterers  but  includes  no 
provision  for  the  interaction  between  scatterers  and  the  Coherent  Potential 
Approximation  (CPA)  which  neglects  the  difference  between  the  external 
excitation  field  and  the  average  field  and  includes  interaction  effects  between 
scatterers.  The  effective  parameters  obtained  using  the  ATA  and  CPA 
approximations  are  contrasted  with  results  obtained  using  several  other  mixing 
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models.  A  detailed  discussion  of  the  approximation  methods  applied  herein  is 
presented  in  Appendix  A. 

4.0  DIELECTRIC  CONSTANTS  OF  FOREST  COMPONENTS 

The  relative  dielectric  constant  effectively  relates  the  polarization  and 
conduction  effects  to  field  quantities  inside  and  outside  an  object.  It  is  usually  a 
complex  quantity  and  can  be  expressed  as, 

f  =  €'  +  jc"  ,  (8) 

where  t'  is  the  relative  permittivity  and  t"  is  the  loss  factor  resulting  from 
propagation  in  a  lossy  medium.  The  loss  factor  accounts  for  losses  arising  from 
polarization  phenomena  in  the  medium  (an  exp(-jwt)  time  dependence  of  the  fields 
is  assumed). 

Application  of  the  approximation  methods  used  in  this  study  to  determine  the 
effective  parameters  of  a  vegetation  environment  requires  a  priori  knowledge  of  the 
dielectric  properties  of  the  individual  scatterers  embedded  within  the  host  medium, 
as  well  as  those  of  the  host  medium  itself  (air  in  this  study).  Experimental 
measurement  of  e  of  the  materials  that  make  *up  a  natural  terrain  or  forest  are 
scarce20,21,  although  several  models  describing  these  values  exist  in  the 
literature2,22,23.  The  values  of  c  are  known  to  vary  with  respect  to  the  moisture 
content,  type  of  material  (i.e.  wood,  leaves,  etc.),  ambient  air  temperature,  wood 
grain  orientation  to  the  electric  field  polarization,  as  well  as  se.  era!  other 
parameters.  Calculation  of  the  effective  parameters  presented  here  are 
representative  of  a  nominal  hard  wood  vegetation  media.  The  relative  dielectric 
constants  of  the  forest  components  were  estimated  using  models  that  calculate 
these  values  as  a  function  of  the  physical  characteristics  of  the  components  (i.e. 
moisture  content,  salinity,  etc.)  and  the  environmental  conditions.  It  was  the  lack 
of  reliable  experimental  data  reporting  measured  values  of  the  dielectric  constants 
of  the  components  of  vegetation  media  that  prompted  the  use  of  models  to 
determine  these  values.  When  experimental  values  of  the  dielectric  constants  of 
foliage  components  for  similar  conditions  to  those  used  in  this  report  were 
available,  they  were  used  for  a  comparison  with  the  model  results.  The  values  of  c 
used  in  this  analysis  correspond  to  values  determined  for  specific  conditions  with 
the  intent  of  showing  common  values  that  can  be  expected  in  a  foliage  medium, 
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and  are  not  professed  to  be  universal  for  all  forest  stands.  The  intention  was  to 
estimate  values  of  e*  for  specific  vegetation  environments  using  nominal  values  for 
the  vegetation  components.  This  would  require  the  reader  to  vary  these  values  in 
order  to  apply  the  formalism  to  a  particular  situation  of  interest.  Detailed 
discussion  of  the  models  used  to  determine  t  for  the  various  vegetation  components, 
and  the  quantitative  results  of  these  calculations,  are  presented  in  Appendix  B. 

5.0  CALCULATION  OF  EFFECTIVE  DIELECTRIC  CONSTANTS 

The  suitability  of  the  approximation  methods  to  estimate  the  effective 
dielectric  constants  of  specific  foliage  media  is  first  examined  by  comparing  values 
of  e*  calculated  using  the  models  with  similar  values  reported  in  the  literature12,24. 
A  computer  program  that  calculates  the  values  of  e*  by  applying  the  mixing 
models  discussed  in  Appendix  A  was  developed  on  a  DEC  VAX9000  computer 
system.  The  program,  in  its  present  form,  produces  plots  of  Re(£*)  and  Im(c*) 
versus  the  volume  fraction  of  an  embedded  scatterer.  If  there  are  multiple 
embedded  scatterers  comprising  the  heterogeneous  medium,  the  model  is  capable  of 
varying  the  volume  fraction  of  only  one  of  the  embedded  scatterers  in  any  one 
calculation.  Because  in  most  cases  wood  will  represent  up  to  95%  of  the  total 
embedded  scatte  :s  volume  fraction  in  a  vegetation  environment,  the  wood 
volume  fraction  was  the  scatterer  varied  in  all  of  the  example  calculations 
presented. 

Model  estimates  of  the  effective  permittivity  and  loss  factor  for  a  forest 
environment  composed  of  hardwood  tree  types  with  the  wood  grain  perpendicular 
to  the  electric  field  polarization  and  ambient  temperature  of  23°C  are  shown  in 
Figures  1  and  2,  respectively.  The  plots  depict  the  approximated  values  of  t*  of  a 
forest  environment  for  an  incident  radiation  frequency  of  1.3GHz.  The  dielectric 
constants  of  the  wood  and  leaves  comprising  the  vegetation  were  obtained  using  the 
values  discussed  in  Appendix  B.  The  volume  fraction  of  the  leaves  in  Figures  1 
and  2  was  assumed  to  be  0.001,  and  the  volume  fraction  of  the  wood  was  varied 
from  0.0001  to  0.1.  The  plot  contains  five  curves  showing  the  results  of  calculations 
using  four  different  shaped  scatterers;  spherical  (ATA  and  CPA),  needles  (Polder), 
cylinders  (Rayleigh),  and  disks  (Bruggeman).  It  is  apparent  that  the  estimations 
provided  by  the  ATA  and  CPA  show  good  agreement  with  one  mother  while 
agreement  between  the  models  defined  for  different  shapes  of  embedded  scatterers 
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Figure  2.  Effective  loss  factor  Im(c*)  of  living  hard  wood  vegetation  media  with  the 
wood  grain  oriented  perpendicular  to  the  electric  field  polarization  Q  23’ C. 
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differ  widely  as  would  be  expected.  The  cylinder  shaped  scatterers  produce  the 
lowest  values  of  c*  while  the  disk  scatterers  produce  the  highest.  The  magnitude  of 
t*  using  needle  shaped  scatterers  is  in  between  the  approximations  for  spherical  and 
disk  shaped  scatterers. 

Model  estimates  of  the  effective  permittivity  and  loss  factor  for  a  forest 
environment  composed  of  hardwood  tree  types  with  the  wood  grain  parallel  to  the 
electric  field  polarization  at  a  temperature  of  23°C  are  shown  in  Figures  3  and  4, 
respectively.  The  plots  depict  the  approximated  values  of  e*  for  the  same 
environmental  conditions  as  presumed  in  figures  1  and  2.  Observation  of  the 
parallel  wood  grain  curves  reveal  that  e*  retains  the;  same  relative  magnitudes  for 
differing  scatterer  shapes  as  that  of  the  perpendicular  wood  grain  curves  with  the 
disk  model  returning  the  highest  estimate  and  the  cylinder  scatterer  model  the 
lowest. 

Comparison  of  the  perpendicular  and  parallel  wood  grain  curves  reveals  that 
the  values  of  Re(e*)  for  the  parallel  wood  grain  tire  consistently  greater  than  those 
of  the  perpendicular  wood  grain.  The  increase  is  only  slight  for  the  sphere  and 
cylinder  curves  and  is  much  greater  for  the  disk  and  needle  curves.  For  the  plots  of 
Im(e*),  the  spherical  and  cylinder  perpendicular  curves  are  consistently  greater 
than  the  corresponding  parallel  curves  versus  volume  fraction,  while  the  difference 
between  the  disk  and  needle  perpendicular  Im(c*)  curves  and  the  corresponding 
parallel  curves  appears  to  increase  with  increasing  volume  fraction,  with  the 
parallel  Im(€*)  values  being  greater.  The  variation  in  the  relative  magnitudes  of  e* 
based  on  electric  field  orientation  are  expected  for  the  different  embedded  scatterer 
shapes.  A  disk  or  needle  would  most  certainly  scatter  differing  amounts  of  energy 
based  on  its  orientation  relative  to  the  electric  field  polarization,  while  spherical 
scatterers  would  scatter  the  same  amount  of  energy  regardless  of  orientation. 

Effective  dielectric  constant  versus  volume  fraction  curves  (not  shown) 
similar  to  those  in  figures  1  through  4  were  generated  for  a  variety  of  incident 
radiation  frequencies  and  several  different  values  of  the  dielectric  constant  of  the 
embedded  scatterers  (representing  different  wood  types,  moisture  content  of  wood, 
salinity,  etc.).  All  of  these  curves  demonstrated  the  same  relationship  between  the 
magnitudes  of  e*  relative  to  the  shapes  of  the  embedded  scatterers  and  wood  grain 
orientation  to  the  electric  field  polarization.  The  changes  in  the  values  of  e* 
relative  to  alteration  of  the  dielectric  constants  of  the  embedded  scatterers  resulted 
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in  only  subtle  changes  in  values  of  Re(c*)  and  somewhat  more  pronounced  changes 
in  Im(c*).  The  dominant  variable  affecting  the  values  of  c*  is  the  volume  fraction 
of  the  embedded  scatterers. 

6.0  COMPARISON  OF  MODEL  PREDICTION  WITH  ESTABLISHED  VALUES 

The  magnitude  of  the  effective  parameters  of  foliage  media  reported  in  the 
literature12,20’25  (where  the  effective  parameters  were  on  the  order  of  1.1+jlO'4) 
correspond  best  with  the  values  of  e*  calculated  using  the  spherical  embedded 
scatterers  (i.e.  the  ATA  and  CPA).  The  values  obtained  by  the  disk  and  needle 
shaped  embedded  scatterer  models  were  much  larger  than  those  generally  reported 
and  ranged  from  1.07  <  Re(t*)  <  2.2  to  1.9(10'2)  <  Im(f*)  <  2.0(10'1)  and 
1.04  <  Re(e*)  <  1.7  to  8.5(10‘3)  <  Im(f*)  <  1.0(10'1)  for  volume  fractions  of 
.001  <  p  <  .05,  respectively.  The  values  obtained  by  the  cylinder  shaped  embedded 
scatterer  models  were  smaller  than  those  reported  in  the  literature  and  ranged  from 
1.006  <  Re(t*)  <  1.1  and  8.1(10'5)  <  Im(f*)  <  1.0(10'3)  for  volume  fractions  of 
.001  <  p  <  .05.  The  major  discrepancy  appears  in  the  complex  term  of  e*  which  can 
vary  by  an  order  of  magnitude  or  more  over  the  range  of  volume  fractions 
examined. 


volume 

fraction 

embedded  scatterer  shape 

sphere 

needle 

cylinder 

disk 

0.001 

1,426 

29 

2,973 

15 

0.01 

293 

14 

569 

7 

0.02 

158 

9 

308 

5 

0.05 

68 

5 

124 

3 

Table  1.  Skin  depth  penetration  in  feet  of  a  l.SGHz  wave  propagating  in  vegetation. 


To  provide  a  tangible  measure  of  the  effect  that  the  variations  in  the 
magnitude  of  e*  has  on  the  approximation  of  the  propagation  loss,  estimates  of  the 
skin  depth  penetration  were  calculated  using  the  different  approximations  of  €*. 
Eq.(5)  was  used  to  determine  the  skin  depth  penetration  of  a  1.3GHz  wave  in 
foliage  media  using  the  values  of  t*  reported  in  figures  1  and  2.  Table  1  presents 
the  results  of  these  calculations  for  various  density  foliage  environments  from  very 
sparse  (0.1%)  to  very  dense  (5.0%).  All  values  shown  in  table  1  reflect  propagation 
through  foliage  when  the  wood  grain  is  perpendicular  to  the  electric  field 
polarization  (i.e.  horizontal  polarization).  Examination  of  table  1  reveals  that  the 
needle  and  disk  models  produce  effective  parameters  that  greatly  under-estimate 
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the  experimentally  measured  values  of  the  skin  depth  penetration  in  foliage  media, 
which  is  reported  to  be  between  50  and  100  feet  in  dense  foliage26,27,  while  the 
spherical  and  cylinder  approximations  provide  a  more  realistic  estimate  of  the 
penetration  depth.  The  skin  depth  penetration  calculated  using  the  effective 
parameters  returned  by  the  CPA  resulted  in  a  penetration  of  300  feet  in  a  normal 
vegetation  medium  and  70  feet  in  a  very  dense  vegetation  medium,  while  the 
effective  parameters  based  on  the  cylinder  model  produced  skin  depth  penetration 
distances  that  were  roughly  twice  those  of  the  spherical  model  estimates. 
Experimental  measurements  of  the  propagation  of  a  1.3GHz  wave  in  foliage 
media28,29  verify  that  the  estimates  of  the  spherical  model  represent  the  best 
approximation  of  the  true  attenuation. 

Agreement  between  the  CPA  and  AT  A  techniques  introduces  the  question 
of  how  important  the  multiple  scattering  processes  are  in  a  vegetation  medium. 
Since  the  AT  A  includes  no  provision  for  the  interaction  between  scatterers  while 
the  CPA  does,  if  the  multiple  scattering  interaction  were  a  dominant  process 
within  the  medium  there  should  be  a  noticeable  difference  between  these 
approximation  methods.  Since  there  is  little  or  no  difference,  it  can  be  concluded, 
based  on  the  results  of  the  model  calculations,  that  multiple  scattering  interactions 
within  the  vegetation  medium  at  certain  frequencies  are  nearly  negligible. 

It  may  also  be  concluded  from  observation  of  figures  1  and  2  that  the  value 
of  e*  is  extremely  dependent  on  the  volume  fraction  of  the  embedded  scatterers 
which  is  representative  of  the  density  of  the  foliage.  In  some  cases  the  magnitude 
of  Im(c*)  can  increase  two  orders  of  magnitude  for  a  sparse  versus  dense  forest 
environment.  The  magnitude  of  Re(e*)  is  also  dependent  on  the  volume  fraction  of 
the  embedded  scatterers  but  to  a  much  lesser  degree  (also  see  Table  1). 

6.1  Comparison  of  Model  Prediction  with  Experimental  Data 

As  another  measure  of  the  validity  of  the  effective  parameters  calculated,  a 
comparison  with  experimentally  measured  values  of  the  attenuation  coefficient  (o 
in  dB/m)  reported  in  the  literature  was  done.  Experimental  measurement  of  the 
attenuation  coefficients  of  EM  waves  propagating  in  vegetation  media  reported  in 
the  literature  are  rare.  Four  sources  were  found  that  contained  measurements  of 
the  attenuation  coefficients  that  were  used  for  comparison  to  the  data  from  this 
analysis. 
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First,  a  comparison  of  the  calculated  attenuation  coefficients  (a  in  dB/m) 
based  on  the  effective  parameters  of  the  four  embedded  scatterer  shape  models  was 
done.  The  results  of  this  calculation  were  plotted  versus  the  frequency  of  the 
incident  radiation  and  are  shown  in  figure  5.  These  values  represent  horizontally 


CALCULATED  ATTENUATION  COEFFICIENTS 


Horizontal  polarization.  1  7.  volume  fraction 


Frequency  (MHz) 


Figure  5.  Relationship  between  a  (dB/m)  calculated  for  various  assumed  embedded 
scatterer  shapes  for  a  1.3GHz  wave  in  a  normal  density  foliage  environment. 


polarized  waves  propagating  through  a  normal  density  foliage  medium  at  am 
ambient  air  temperature  of  23°C.  As  with  the  skin  depth  calculation,  the  models 
based  on  the  disk  and  needle  embedded  scatterers  produced  attenuation  coefficients 
much  greater  than  those  obtained  by  the  spherical  and  cylinder  based  models. 
Determination  of  which  of  these  models  provides  the  best  estimate  of  the 
propagation  loss  is  verified  by  comparison  of  the  cadculated  attenuation  coefficients 
with  experimentally  measured  vadues.  From  this  comparison  it  was  concluded  that 
the  effective  parameters  calculated  using  the  CPA  provided  the  best  estimate  of  o. 
Figures  6  and  7  present  a  comparison  of  the  meatsured  amd  cadculated  attenuation 
coefficients  for  horizontal  and  vertical  polarized  waves,  respectively. 

Figure  6  depicts  the  attenuation  coefficients  calculated  using  the  CPA  for 
embedded  scatterer  volume  fractions  of  1%  and  5%  compared  to  experimental 
measurements  of  these  values  from  four  sources28,29,30.  The  measured  vadues 
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Figure  6.  Comparison  of  model  estimates  of  a  (dB/m)  with  experimentally 
measured  values  for  dense  foliage.  Horizontal  polarization. 


Figure  7.  Comparison  of  model  estimates  of  a  (dB/m)  with  experimentally 


measured  values  for  dense  foliage.  Vertical  polarization. 


represent  the  attenuation  loss  of  horizontally  polarized  waves  propagating  in  dense 
foliage.  The  calculated  values  appear  to  provide  a  good  estimate  of  the  attenuation 
as  the  majority  of  the  experimentally  measured  points  fall  between  the  1%  and  5% 
curves.  As  would  be  expected  the  experimental  values  are  near  the  upper  curve 
because  all  of  the  measured  values  represent  propagation  through  what  was 
described  as  very  dense  foliage.  It  is  apparent  from  the  figure  that  the  best 
agreement  between  the  model  estimates  and  experimental  measurements  occurs 
between  the  frequency  limits  of  50MHZ  to  1GHz.  In  this  range  the  experimentally 
measured  values  appear  to  fall  within  the  1%  to  3%  volume  fraction  range  expected 
of  dense  foliage.  At  lower  frequencies  the  model  appears  to  agree  with  some 
measurements  while  underestimating  others.  At  frequencies  above  1GHz,  the 
model  appears  to  provide  reasonable  estimates  of  the  attenuation  even  though  the 
frequency  limitations  of  the  mixing  models  begin  to  become  apparent  somewhere  in 
the  1GHz  to  3GHz  range.  There  is  a  large  disparity  between  experimentally 
measured  values  of  the  attenuation  coefficients  at  certain  frequencies.  In  some 
cases  this  disparity  is  greater  than  the  difference  between  the  model  prediction  and 
measured  data.  The  difference  between  the  experimentally  measured  and  model 
predicted  values  could  be  attributed  to  any  one  of  several  possible  factors  that 
could  contribute  to  natural  variation  in  the  dielectric  properties  of  the  vegetation 
media  in  which  the  measurements  were  taken.  Essentially,  variation  in  the 
dielectric  properties  of  the  components  comprising  different  vegetation 
environments,  and  even  within  vegetation  environments,  prevent  consistently 
accurate  determination  of  the  attenuation  properties  of  vegetation  media  by 
experimental  or  theoretical  means. 

Figure  7  shows  the  attenuation  coefficients  calculated  using  the  CPA  for 
embedded  scatterer  volume  fractions  of  1%  and  5%  compared  to  experimented 
measurements  of  these  values  by  three  sources.  The  measured  values  represent  the 
attenuation  loss  of  vertically  polarized  waves  propagating  in  very  dense  foliage. 
The  values  calculated  based  on  the  model  do  not  provide  as  good  an  estimate  of  the 
attenuation  of  the  vertically  polarized  waves  as  for  the  horizontally  polarized 
waves.  Although  several  of  the  experimented  points  fall  on  or  near  the  5%  curve, 
the  majority  of  the  experimentally  measured  points  are  above  the  5%  curves.  True 
volume  fraction  values  would  not  be  expected  to  be  much  above  1%  to  3%  in 
practical  terms,  which  means  that  the  model  will  underestimate  the  attenuation 
loss  for  vertically  polarized  waves.  The  greater  underestimation  of  the  attenuation 


15 


coefficients  of  the  vertically  polarized  waves  than  the  horizontally  polarized  waves 
can  be  attributed  to  the  greater  depolarization  effects  that  vertical  waves 
experience  compared  to  horizontal  waves  when  propagating  through  vegetation30. 
This  reveals  an  inherent  weakness  in  the  effective  medium  approximation 
technique  which  does  not  account  for  polarization  related  scattering  effects 
resulting  from  differently  shaped  embedded  scatterers  in  the  foliage  and  relies 
solely  on  the  variation  of  the  dielectric  constants  of  the  embedded  scatterers  to 
account  for  polarization  effects. 

The  general  shapes  of  the  curves  in  figures  6  and  7  appear  to  mimic  the 
trends  exhibited  by  the  reported  experimental  data  versus  frequency.  The 
magnitude  errors  could  be  attributed  to  poor  estimates  of  the  dielectric  constants  of 
the  embedded  scatterers  comprising  the  foliage  medium  compared  to  the  true 
values  of  the  forests  where  the  propagation  measurement  were  taken.  Dielectric 
constants  of  the  foliage  components  can  be  expected  to  vary  from  forest  to  forest  as 
well  as  spatially  within  a  forest.  Therefore,  better  estimates  of  the  composition  of  a 
vegetation  media  will  lead  to  better  approximations  of  the  dielectric  constants  of 
the  components  which  will  result  in  better  estimates  of  the  effective  parameters. 

Overall,  the  effective  dielectric  constants  calculated  using  the  CPA  will 
provide  a  good  estimate  for  of  the  attenuation  loss  of  a  horizontally  polarized  wave 
and  a  fair  estimate  of  a  vertically  polarized  wave  propagating  in  a  vegetation 
media.  The  frequency  bounds  are  difficult  to  verify  based  on  the  results  shown  in 
figures  6  and  7.  Given  the  accuracy  of  the  CPA  model  estimates  of  the  effective 
parameters,  use  of  the  other  approximation  techniques  would  clearly  result  in 
either  much  to  great  an  attenuation  loss,  in  the  case  of  the  needle  and  disk  models, 
or  much  to  little  attenuation,  in  the  case  of  the  cylinder  based  model. 

6.2  Application  of  the  Dielectric  Slab  Model 

As  a  final  test  of  the  suitability  of  the  effective  parameters,  calculated  using 
the  CPA,  to  model  attenuation  through  vegetation,  the  results  of  the  application  of 
the  slab  model31  in  a  scattering  simulation  are  shown.  In  this  simulation,  a 
transmitting  antenna  and  receiving  antenna  are  located  over  a  surface  with  two 
scales  of  roughness  that  is  covered  by  a  foliage  layer.  Both  the  coherent  and 
incoherent  power  scattered  into  the  receiver  tire  calculated.  The  model  assumes 
only  in-plane  scattering;  there  is  no  out-of-plane  aziumthally  dependent  scattering. 
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Waves  are  launched  from  the  transmitting  antenna  down  to  the  foliage  covered 
terrain.  Each  ray  from  the  transmitting  antenna  is  treated  as  a  plane  wave  whose 
amplitude  is  weighted  by  the  pattern  factor  of  the  transmitting  antenna.  Each 
wave  in  the  angular  spectrum  of  plane  waves  emanating  from  the  transmitter  is  ray 
traced  through  the  foliage  layer  down  to  the  rough  surface,  is  reflected,  and  then 
ray  traced  back  up  through  the  foliage  layer  and  into  a  bistatic  receiving  antenna. 
The  rough  surface  scattering  is  governed  by  the  two  scales  of  roughness  model  of 
Peake-Barrick31.  Only  the  mean  value  of  the  electric  field  (coherent  part)  is  ray 
traced  down  through  the  foliage  layer.  The  coherent  wave  reflected  from  the  rough 
surface  is  ray  traced  back  up  through  the  foliage  layer  and  is  amplitude  weighted 
by  the  pattern  factor  of  the  receiving  antenna.  Also,  the  diffuse  EM  energy 
reflected  from  the  rough  surface  is  composed  of  an  angular  spectrum  of  waves,  each 
of  which  is  ray  traced  up  through  the  foliage  layer,  suffers  some  attenuation,  and  is 
then  amplitude  weighted  by  the  pattern  factor  of  the  receiving  antenna.  In  this 
investigation,  the  theoretical  model,  with  appropriate  rough  surface  parameters, 
was  used  to  predict  an  effective  normalized  cross  section  a  which  is  then  compared 
with  experimental  data.  The  comparisons  shown  here  are  for  three  terrain  types:  a 
grass  covered  field,  a  forest  and  a  flooded  forest.  The  agreement  between  the 
theoretical  model  and  the  data  is,  in  general,  quite  good. 

Figures  8  and  9  show  L-Band  scattering  data  (cr°  vs.  0;)  for  a  1-meter  grass 
field  taken  with  a  SAR32.  The  theoretical  model  calculations  were  performed  with 
the  following  surface  parameters:  <rL=1.22  m,  TL=5  m,  <r,=3.16  •  10' 2  m, 
T*=4.0  •  10'2m,  where  c(ground)  =  30.0+j0.6  and  e*(grass)=1.2+j0.005.  Figure  8 
shows  the  excellent  agreement  between  the  SAR  data  and  the  theoretical  model  for 
HH  polarization  for  the  grass  covered  field.  Figure  9  shows  moderately  good 
agreeme'  ‘  between  theory  and  data  for  the  VV  polarization  for  the  grass  covered 
field.  The  5  dB  discrepancy  may  be  due  to  some  depolarization  for  the  VV 
polarization  (the  vertically  polarized  incident  EM  waves  axe  partly  scattered  into 
horizontally  polarized  waves). 

Figure  10  shows  two  data  points  (cr°)  for  a  sparse  forest  (  ~  22  m  trees)  from 
a  SAR  with  HH  polarization.  One  data  point  (the  solid  circle)  corresponds  to  a 
normal  forest32,  and  falls  very  near  the  theoretical  model  curve.  The  surface 
parameters  used  in  the  theoretical  model  curve  are  given  as  follows:  <rL=  1.22m, 
TL=10.0m,  <r,=3.16  •  10’2m,  and  Tg=4.0  •  10'2m,  while  t(ground)=  30.0+j0.6  and 
e*(forest)  =  1.03+j0.0002.  The  second  data  point  (the  empty  circle)  in  Fig.  10 
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ONE  METER  GRASS  COVERED  FIELD 
L-Band,  horizontal  polarization 


Figure  8.  Comparison  of  model  predicted  a  *  values  with  measured  SAR  data  for  a 
1  meter  grass  field  at  L-Band  and  horizontal  polarization. 


Figure  9.  Comparison  of  model  predicted  a  °  values  with  measured  SAR  data  for  a 
1  meter  grass  field  at  L-Band  and  vertical  polarization. 
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Figure  10.  Comparison  of  model  predicted  a  *  values  with  measured  SAR  data  for  a 
22  meter  forest  at  L-Band  and  horizontal  polarization. 


corresponds  to  a  for  the  same  forest  (  ~  22  m  trees)  flooded.  The  theoretical  model 
surface  parameters  are  the  same  as  for  the  normal  forest,  but  with  t(ground)= 
80.0-fj0.6  (corresponding  to  lake  water).  The  agreement  of  the  data16  with  the 
theoretical  model  is  quite  good  (within  3dB).  The  theoretical  model  also  shows  the 
same  trend  as  the  data  points,  the  <r°  for  the  flooded  forest  is  several  dB  higher  than 
a  for  a  normal  forest. 

In  Fig.  11,  several  L-Band  data  points  for  a  from  a  forest  comprised  of  10m 
trees  are  shown  versus  incidence  angle  0;  for  HH  polarization33.  The  surface 
parameters  used  in  the  theoretical  calculations  are  as  follows:  oL=  1.22m, 
TL=10.0m,  <rs=3.16-  10'2m,  T*=4.0- 10"2m,  where  c(ground)=30.0+j0.6  and 

c(forest)=1.03+j0.0002.  The  agreement  of  the  data  with  the  theoretical  model  is 
again  quite  good.  The  most  statistically  significant  point  of  the  measured  data 
corresponds  to  the  value  of  <r°  at  an  incidence  angle  of  55’.  This  is  also  where  the 
theoretical  model  curve  nearly  intersects  the  data  curve. 
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10  METER  SPARSE  FOREST 

L-Bond,  horizontal  polorization 


Figure  11.  Comparison  of  model  predicted  <r°  values  with  measured  data  for  a  10 
meter  forest  at  L-Band  and  horizontal  polarization. 

7.0  CONCLUSION 

Based  on  the  overall  agreement  between  the  model  approximations  and 
experimental  data,  the  effective  dielectric  constants  calculated  using  the  Average 
T-Matrix  Approximation  and  the  Coherent  Potential  Approximation  provide  the 
basis  for  a  reasonable  formalism  to  predict  the  propagation  loss  of  an 
electromagnetic  wave  propagating  through  a  vegetation  environment.  Values  of 
the  effective  dielectric  constants  calculated  agree  well  with  the  values  currently- 
accepted  in  the  literature.  In  addition,  the  calculated  attenuation  coefficients  of 
waves  propagating  through  vegetation  media  and  experimental  measurements  of 
these  values  agree  quite  well.  Finally,  as  demonstrated  by  the  bistatic  scattering 
model  output  (figures  8-11),  the  dielectric  slab  model  will  provide  a  good  estimate 
of  the  amplitude  loss  of  a  coherent  electromagnetic  wave  propagating  through 
foliage.  The  formalism  is  not  computationally  intensive  and  yet  accounts  for  the 
effects  of  multiple  scattering  phenomenon  and  is,  therefore,  appealing  for  use  in 
large  scale  system  simulation  programs  where  computation  time  has  to  be 
minimized.  Use  of  this  formalism  as  the  primary  method  to  determine  the 
propagation  effects  of  a  wave  propagating  great  distances  through  foliage  media  is 
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not  recommended.  The  model  is  primarily  intended  to  provide  a  secondary 
attenuation  factor  attributable  to  the  presents  of  vegetation.  That  is,  to  be  applied 
in  a  system  simulation  p.ogram  similar  to  the  bistatic  scattering  model  discussed  in 
this  report. 

The  primary  variable  influencing  the  magnitude  of  the  effective  par  angers 
is  the  volume  fraction  of  the  embedded  scatterers  which  is  proportional  to  the 
density  of  the  foliage.  This  effect  appears  to  be  consistent  with  the  experimental 
measurements  reported  in  the  literature  which  do  vary  based  on  the  type  and 
condition  of  the  vegetation  but  are  mainly  dependent  on  the  density  of  the  foliage. 

Because  of  inherent  limitations  in  the  mix' \g  model,  the  size  of  the 
embedded  scatterers  as  compared  to  Jhe  i  ncident  radiation  frequency  will  influence 
the  accuracy  of  the  model  estimates  of  propagation  loss.  Therefore,  application  of 
the  effective  parameters  will  1  ,  best  suited  to  model  propagation  through  a  forest 
canopy  or  vegetation  having  snail  diameter  trunks,  each  of  which  present  smaller 
dimension  scatterers  than,  say,  the  trunks  of  a  mature  forest.  Another  important 
effect  of  the  wavelength  to  embedded  scatterer  size  is  the  frequency  limitations  of 
the  model.  It  appears  from  the  data  that  the  model  will  provide  good  estimate  of 
the  attenuation  loss  up  to  frequencies  of  3.2GHz  in  certain  vegetation 
environments. 

As  shown  in  figures  6  and  7,  the  ability  to  account  for  specific  polarization  is 
weak.  This  can  be  primarily  attributed  to  the  lack  of  variation  in  the  scattering 
properties  of  the  embedded  scatterers,  which  are  assumed  to  be  constant  geometric 
shapes.  The  variation  in  polarization  that  is  demonstrated  is  entirely  attributable 
to  the  differing  dielectric  constants  of  the  embedded  scatterers  based  on  the 
orientation  of  the  wood  grain  to  the  electric  field  polarization.  In  addition  to  the 
inability  to  accurately  distinguish  the  attenuation  for  differently  polarized  waves, 
the  proposed  formalism  is  only  capable  of  determining  the  amplitude  loss  of  the 
wave  and  includes  no  provisions  to  determine  the  phase  of  the  wave.  If  necessary, 
the  model  could  be  expanded  to  include  the  effects  of  different  scattering  shapes 
and  could  possibly  provide  a  better  approximation  of  the  attenuation  of  waves 
having  different  polarizations.  At  present,  the  approximation  of  the  attenuation 
loss  is  sufficient  for  the  applications  for  which  this  formalism  was  developed  and 
the  possible  improvements  gained  by  the  expansion  of  the  model  is  not  considered 
a  necessity. 
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APPENDIX  A 


MODELS  FOR  THE  EFFECTIVE  DIELECTRIC  CONSTANTS  OF  FOLIAGE 

In  this  section  a  review  of  the  expressions  used  to  calculate  the  effective 
dielectric  constants  e*  of  foliage  media  is  presented.  Because  the  derivation  of 
these  equations  are  all  readily  available  in  the  literature5,9,16,18  the  derivations 
shown  here  will  be  kept  to  a  minimum  in  the  interest  of  space. 

Subject  to  the  constraints  stated  in  the  body  of  the  paper,  the  continuous 
medium  approach,  characterized  by  an  effective  dielectric  constant,  may  be 
analyzed  by  solving  the  wave  equation  with  a  fluctuating  (random)  dielectric 
constant  using  various  approximation  techniques.  The  transition  operator  T  is 
used  as  a  basis  for  the  various  approximations  in  multiple  scattering  formalisms. 
Maxwell’s  equations  are  written  as, 


and, 


V  •(<„!)  +  £  (<.-<<,)  2  v.(*E)=o, 

1=1  J=1  1 

VxE  =  0 


(1) 

(2) 


where  (2)  is  a  quasi-static  approximation.  Use  of  the  quasi-static  approximation 
assumes  that  the  dimensions  of  the  scatterers  are  much  less  than  the  wavelength  of 
the  incident  radiation  thus  imposing  an  upper  frequency  limitation  on  the  results  of 
the  approximation  models.  The  frequency  limitations  are  governed  by  Eq.  (4)  in 
the  main  body  of  the  report,  and  axe  dependent  on  the  frequency  of  the  incident 
radiation,  effective  dielectric  constant  of  the  medium  and  dimension  of  the 
scatterers.  Equation  1  implies  that  there  axe  m  constituents  and  N-  is  the  number 
of  scatterers  of  the  i^  type.  Here,  pi  is  the  volume  fraction  of  the  i^  type  of 
scatterer,  x  is  a  spatial  parameter  and  is  the  support  of  the  scatterer  of  type 
i,  where, 


1  x  in  jtl1  scatterer  of  type  i 
0  otherwise 


(3) 


The  two  relationships, 


D(x)  =  e(x)E(x) 

and, 


(4) 
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define  c*  where  <  •  >  denotes  ensemble  average. 
Expressing  Eq.  (1)  in  operator  form, 


(Lo+M)E  =  0  ,  V  xE  =  0  ,  (6) 

where, 

Lo=V -(<„•),  (7) 

M  =  M|+  ...  +  Mm, 

Wlth7  N. 

M,  =  E  V0  ,  (8) 

and, 

V.j=  («r  «o)’vUy(*))  •  (9) 


The  T  matrix  may  be  defined  as, 

T  =  (L„  +  M)_1M  .  (10) 

After  considerable  algebraic  manipulation,  and  expanding  the  T  matrix  in  a  power 
series  and  neglecting  higher  order  multiple  interaction  terms,  an  expression  for  <* 
may  be  determined  in  terms  of  expectation  values  of  the  T  operator  and  the 
dielectric  constant  ’  . 


~  [<c>  -  £<eT,>]-[i-  £<t i>yl  (ii) 

1=1  1=1 

In  addition,  if  it  is  assumed  that  there  is  only  a  single  class  of  embedded  scatterers 
(m=l)  and  the  scatterers  axe  spherical  in  shape,  the  following  approximation  for  e* 
may  be  obtained5,9: 


(12) 


where  e0  represents  the  dielectric  constant  of  the  host  medium,  px  and  ej  the 
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volume  fraction  and  dielectric  constant  of  the  embedded  scatterer,  respectively,  and 


<‘>  -  (1  •  PX)‘<,  +  ?!<1 


(13) 


Spherical  embedded  scatterers,  although  not  representative  of  the  actual  shapes  of 
vegetation  components,  were  used  to  demonstrate  the  effectiveness  of  the  dielectric 
mixing  models  at  estimating  the  values  of  e*.  The  best  physical  analogy  between 
the  spherical  scatterer  and  actual  vegetation  components  would  be  clumps  of  leaves 
and  branches  in  a  vegetation  canopy  that  can  be  somewhat  spherical  or  infinitely 
small  spheres  representing  scattering  centers  at  the  actual  point  of  scattering. 

For  low  volume  fractions,  the  average  T-matrix  approximation  (ATA)  may 
be  written, 


’■('o'-’vf,)  -  +  3<-,(^+^)  +  3 ) 


(14) 


This  approximation  is  good  only  for  small  volume  fractions  p  <  0.02  and  neglects 
interparticle  scattering  interactions. 

In  the  coherent  potential  approximation  (CPA),  the  difference  between  the 
field  exciting  the  medium  and  the  average  field  is  neglected5,9,18,34.  A  reference 
dielectric  constant  cr  is  introduced  into  Maxwell’s  equations,  and  is  later  chosen  to 
simplify  (optimize)  the  equations  for  e*: 


and, 


— >  m  /  — K 

V-(crE)+  E  EV-((erc0)x..Ej  + 
i=l  j=i  '  lj  ' 

VxE  =  0 


(15) 


In  operator  notation,  this  is  written  as, 


where, 


(Lr  +  Mr)E  =  0  , 

(16) 

Lr=V  •  (er  • )  , 

(H) 

m  N  < 

Mr  =  E  EVf,  , 

•=i  j=i 

(18) 
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and, 


v;,=v.((<,-  <»)■*,)  + (-^> 


The  T  matrix  in  this  approximation  is  defined  as, 

€*  ~  [<e>  -  £  E«T,>].[I-  E  E<T,>]1 

»=i  i=i  i=x  i= 1 


where, 

T fj  =  (Lr  +  VfjrWfj 


(19) 


(20) 


(21) 


After  much  algebraic  manipulation,  an  expression  for  e*  may  be  determined  in 
terms  of  er,  and  expectation  values  of  the  dielectric  constant  and  the  Tr  matrix. 
This  expression  for  c*  greatly  simplifies  if  er  is  chosen  so  that5’9’18’34, 


N, 


<  E  LTfj>  =  0 

»= i  j= i 


This  is  the  CPA  and  requires  that, 


£o  -  fr  _  ^  (  <o-tj  \ 

3«r  1=1  •  \3er  +  €--e0J 


(22) 


(23) 


where  p.=^rS3Ci  with  c,  being  the  average  number  of  scattering  centers  of  class  *  per 
unit  volume.  The  general  expression  for  e*  obtained  by  solving  the  operator 
equation  (20)  is, 


e-  ~  e. 


where, 


o+  ^E  Pfr  -  P€oj  -  (!-p)co(%r)  + 


(24) 


p 


When  only  one  class  of  embedded  scatterer  is  considered  (m=l)  and  the  scatterers 
are  very  small,  Eqs.  (23)  and  (24)  combine  to  give  the  CPA  for  e*; 


eo;fi^i)-  €o^l  +  : 


(  £r£o  ^  , 

2/(«r«o)2(f?  +  13£o«r5£o  )\\ 

l\2f0  +  cj  +  P 

*  £o(£i +  2£o)3  7 

(25) 
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Given  the  questionable  correlation  between  the  actual  physical  makeup  of 
vegetation  components  and  spherical  embedded  scatterers,  mixing  formula  derived 
using  various  embedded  scatterer  shapes  (needles,  cylinders,  and  disks)  were  also 
examined.  Following  the  derivation  of  Polder  and  Van  Santen  VanBeek16  the 
generalized  expression  for  the  effective  parameter  of  a  heterogeneous  dielectric 
mixture  is  expressed, 


e-cr 


_  1 


3Vl({rfo).£[^4-o)] 


(26) 


where  a,b,c  represent  the  axes  of  the  embedded  scatterer,  Aj  is  the  depolarization 
factor  along  the  i^“  axis  of  the  embedded  scatterer,  and  e0  is  specified  later  and 
must  lie  within  the  range  of  €0  <  £0  <  f  • 

For  long  prolate  spheroids  A;  =  (5-6,  i-6,  2 S),  where  6  <  1.  Inserting  these 
values  for  Aj  into  Eq.(26)  leads  to, 


e 


1 

€0  +  2£(e1-£o) 


(27) 


Equation  (27)  can  be  solved  for  various  scatterer  shapes.  The  special  case  of  6-»0 
represents  embedded  scatterers  of  long  thin  needles.  Solved  in  the  limits  t0=  c0, 


e 


.  ..  (ex-fo)(5f0  +  €i) 

tl+  1  3(<0-tl) 


(28) 


Another  formula  obtained  by  Rayleigh16  represents  embedded  scatterers 
that  axe  rectangular  cylinders  with  their  axes  parallel  and  perpendicular  to  the  field 
direction.  The  resulting  formula  is  expressed, 

fo  +  crvi(greo) 

°  eo  +  ei  +  vi(€i'€o) 


Disks  can  be  approximately  regarded  as  oblate  spheroids.  For  flat  extended  oblate 
spheroids  Aj  =  (6,  6,  1-2$),  where  6  <  1.  Equation  26  cam  be  applied  again  to  obtain 
a  general  formula  to  estimate  c  for  disk  shaped  scatterers. 


(  —  tr 


t(£r£o)£o  ( 


+  *(£r«o)  £r2$(<rc0) 


) 


(30) 
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For  the  case  when  6=0  and  A;  =  (o,  o,  1)  Eq.(30)  becomes, 


for  e0=  £0.  Bruggeman  was  the  first  to  obtain  Eq.(31)  by  means  of  his  integration 
scheme  which  resulted  in  the  following  equation  for  disks16, 

_  3c0  +  2v1(f1-€0) 

It  is  possible  to  extend  the  case  of  m=l  (one  embedded  scatterer  in  a  host 
medium)  to  the  case  of  severed  embedded  scatterers  in  a  host  medium  (m  >  2)  by 
using  an  iterative  procedure.  When  a  background  medium  £0  contains  the  classes 
of  embedded  scatterers,  and  t2,  the  effective  dielectric  constant  of  the  medium 
can  be  defined  as  e,(c0;e1,p1/(i  -  p2))  where  the  embedded  particles  are  characterized 
by  (c2,^>2).  An  alternate  to  this  would  be  to  determine  the  effective  dielectric  as 
£*(£o5£i ,pj(x  *  Px))  where  the  embedded  scatterers  are  characterized  by  (cj,p  ). 
This  recursive  technique  constructs  models  of  media  containing  several  classes  of 
scatterers  by  embedding  the  "next”  class  of  scatterers  into  a  medium  whose 
effective  parameters  have  been  derived  for  the  "previous”  classes  of  scatterers2. 

For  the  case  of  two  embedded  scatterers  where  the  volume  fractions,  p 1  and 
P2,  are  in  the  low  to  moderate  range, 


£*(£o;£i>Pi;£2>p2)  -  £*(£*(£o;£i^1  'p  ^);£2^2) 

(33) 

£*(£0;£1>P!;£2^2)  »  £*(£*(£0)£2^1  2p  j)£l ,PX) 

(34) 

will  yield  good  approximations  of  e*  of  the  composite  medium.  The  terms  pj{ i  - 
P2)  and  pj{ l  -  pj  in  Eqs.  (33  and  34)  represent  the  effective  volume  fra-  tions  of 
the  intermediate  mixtures  of  ^fj)  and  (f0,c2),  respectively.  When  pt  <  p2,  Eq.  (33) 
will  be  the  preferred  approximation  and  if  p2<  pv  Eq.  (34)  will  be  the  preferred 
approximation. 

To  compute  e*  of  a  medium  with  M  classes  of  embedded  scatterers,  a 
technique  similar  to  "iterated  homogenization”  is  applied35.  One  can  determine  the 
approximation  of  «-*(c0;e1,/>1;  •  .  ;e»n,pm),  where  p.  <  pi+v  by  recursively  adding 
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scatterers  to  the  "new"  background  medium  in  order  of  increasing  density.  The 
effective  volume  fraction  of  a  scatterer  of  type  j  in  the  medium  defined  by  the 
background  medium  and  the  first  j  scatterers  (e0,elt...,c  -)  may  be  expressed, 

0^uk-JU—  .  (35) 


1-Ep, 
1=1  3 


Using  this  expression,  the  estimate  of  em  can  be  expressed  as5, 


£*(£0>£1>Pi>-'>£M’Pm)  ~  £*(£*{£o;cl^1  _  p  .  p  ^)>£M’Pm^ 

~  .  .e*(€*(e0;£x,Pi-M);£2,p2  M}‘  •  ^M-hP^-J^M’Pm)  •  (36) 
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APPENDIX  B 


DIELECTRIC  PROPERTIES  OF  TREES 

A  limited  number  of  studies  have  been  done  that  experimentally  or 
theoretically  determine  the  dielectric  properties  of  trees1’2,3,20,21,23.  Despite  this, 
reliable  dielectric  measurements  of  living  wood  axe  scarce,  although  some  data  has 
been  compiled  for  dry  wood  or  lumber.  The  major  difference  between  the  dielectric 
properties  of  lumber  and  those  of  living  wood  are  attributable  to  the  moisture 
content  of  each.  Some  differences  could  be  due  to  sap  and  other  nutrients  within 
the  wood  that  may  also  affect  the  electrolytic  nature  of  the  wood,  but  these  effects 
are  secondary  when  compared  to  the  moisture  content22. 

The  model  used  in  this  analysis  to  determine  t'  of  wood  was  developed  by 
James11  and  further  enhanced  by  Brown22  and  includes  the  effects  of  incident 
polarization  and  frequency.  Values  of  c'  of  the  hard  and  soft  wood  trees  calculated 
using  the  model  do  not  vary  a  great  deal  over  the  frequency  range  considered  in 
this  report.  Data  reported  from  experimental  measurements  of  t'  of  living  wood 
also  exhibited  a  leveling  off  above  20MHz2,20,21.  These  data  do.  however  show  a 
slight  but  constant  downward  trend  in  the  magnitude  of  c#  as  the  frequency 
increases  from  audio  to  radio  to  infra-red  frequency  bands.  By  close  examination  of 
these  data  it  was  determined  that  t'  consistently  decreased  over  the  frequency 
range  of  interest  in  this  study.  The  values  of  «'  reported  in  Figure  IB  reflect  this 
slight  downward  trend  as  the  frequency  increases.  The  magnitude  of  the  changes 
were  based  on  the  results  of  experimental  data  reported  for  wood  of  similar 
densities  and  moisture  contents20.  The  slight  decrease  in  c'  as  the  frequency 
increases  may  be  attributable  to  the  fact  that  at  higher  frequencies  the  time 
required  for  polarization  to  occur  may  be  equal  to  or  even  exceed  the  period  of 
alternation  of  the  applied  electric  field.  Therefore,  the  magnitude  of  the 
polarization  decreases  resulting  in  a  subsequent  decrease  in  t'. 

The  magnitude  of  e'  of  living  wood  can  be  expected  to  increase  as  the 
ambient  temperature  approaches  the  freezing  point3.  Experimental  data  reported 
seems  to  substantiate  the  slight  increase  in  e'  as  temperature  decreases  but  this 
conclusion  was  not  apparent  from  the  data  in  all  reports  cited  Skaar  49,  Yavorsky 
51,  Brown  82,  and  Broadhurst  70]. 
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Relative  Permittivity  of  Living  Hardwood 
Ambient  Air  Temperature  of  23*  C 


Figure  IB.  Re(c*)  of  living  hard  wood  with  the  incident  radiation  electric  field 
polarization  either  parallel  oriented  or  perpendicular  to  the  wood  grain  orientation. 


LOSS  FACTOR  OF  WOOD 

The  dependence  of  the  loss  factor  ( e ")  on  the  physical  properties  of  wood  is 
more  difficult  to  determine  than  that  of  e'.  A  great  disparity  is  reported  to  exist 
between  the  loss  factor  of  cut  wood  and  a  sample  from  a  living  tree  branch  ’  . 
The  loss  due  to  dipolar  depolarization  is  sometimes  the  only  type  of  loss  considered 
in  high  frequency  models  for  leaves  and  wood.  In  wood  and  leaves,  water  volume  is 
considered  to  dominate  the  variation  of  e"  and  so  a  mixing  type  of  formula  is 
generally  based  on  the  estimates  of  the  water  volume  in  a  unit  volume  of  material. 
A  model  to  determine  the  wood  loss  factor,  accounting  for  each  type  of  dominant 
loss  phenomenon,  has  been  derived  by  Brown22  and  is  frequency-dependent  on  the 
incident  radiation,  the  relaxation  frequency  of  water  and  the  moisture  content  of 
the  wood.  The  results  from  this  model  compared  favorably  with  experimental 
data20,  although  widespread  application  is  doubtful.  The  frequency  dependence  of 
t"  on  the  radiation  is  quite  complicated  and  may  be  considerable  in  magnitude. 
The  variation  of  t"  with  respect  to  the  temperature  has  been  shown 
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experimentally2  to  decrease  in  moist  wood  as  the  temperature  increases.  This 
reversal  is  attributable  to  the  fact  that  at  higher  temperatures  the  water  dipoles 
experience  less  frictional  effects.  The  correlation  of  c"  with  the  wood  grain 
orientation  relative  to  the  electric  field  has  not  been  established  experimentally. 
The  loss  factor  as  measured  with  an  electric  field  parallel  to  the  wood  grain  is 
greater  than  that  with  the  electric  field  perpendicular  to  the  grain  as  the  model 
predicts.  The  reason  for  this  behavior  may  be  attributed  to  the  basic  structure  of 
cellulose  crystallites2. 

Values  of  e"  of  wood  for  conditions  of  interest  in  this  study  axe  shown  in 
Figure  2B  as  a  function  of  frequency.  The  calculations  of  these  hardwood  values 
were  based  on  a  mean  moisture  content  of  80%,  and  wood  density  of  0.75.  It  is 
common  for  the  wood  moisture  content  when  taken  as  a  percentage  of  the  oven 
dried  weight  of  the  wood,  as  defined  by  wood  technologists  and  forestry 
professionals,  to  have  values  greater  than  100%36.  All  values  were  determined 
assuming  an  ambient  air  temperature  of  23°C  ’  . 


Relative  Loss  Factor  of  Living  Hardwood 
Ambient  Air  Temperature  of  23*  C 


Frequency  (Hz) 


Figure  2B.  Im(i*)  of  living  hard  wood  with  the  incident  radiation  electric  field 
polarization  either  parallel  oriented  or  perpendicular  to  the  wood  grain  orientation. 
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COMPLEX  RELATIVE  DIELECTRIC  CONSTANTS  OF  LEAVES 


To  obtain  a  valid  estimate  of  the  effective  parameters  of  a  foliage  medium, 
values  of  the  dielectric  properties  of  leaves  must  be  known.  Although  some 
reports3  site  that  the  leaves  on  leaf-bearing  trees  can  be  ignored  below  frequencies 
of  3GHz,  the  scope  of  this  study  includes  operating  frequencies  near  this  boundary 
which  will  necessitate  the  inclusion  of  the  effects  of  leaves.  There  are  very  few 
sources  in  the  literature  that  provide  measured  data  for  e  of  leaves.  Those  that  are 
available20,21  report  e  of  only  certain  types  of  leaves  and  for  a  limited  number  of 
incident  radiation  frequencies.  To  obtain  approximations  for  these  values,  a 
modeling  technique  developed  by  DeLoor37  that  determines  e  of  leaves  based  on 
their  moisture  content,  the  ambient  air  temperature,  and  the  incident  radiation 
frequency  was  used.  In  the  DeLoor  model  the  leaf  is  modeled  by  considering  the 
leaf  to  be  comprised  of  a  background  medium  composed  of  organic  material,  and 
that  ellipsoidal  disks  of  water  are  randomly  distributed  throughout  this  material. 
The  model  uses  a  dielectric  mixing  formula  that  is  a  function  of  the  volume 
fraction  of  the  water  inclusions  and  e  of  the  host  material  and  water. 

The  relative  dielectric  constants  of  leaves  for  various  conditions  of  interest 
in  this  study  have  been  computed  and  are  shown  in  Table  IB.  Assuming  that  the 
moisture  content  of  the  leaves  is  high  for  the  tropical  environment  in  which  much 
of  the  experimental  measurements  of  propagation  through  foliage  were  done,  a 
moisture  content  of  65%  is  used  in  the  determination  of  the  dielectric  constant  of 
leaves  in  the  study.  The  salinity  of  leaves  can  vary  from  1%  to  12%.  As  there  is 
no  way  to  estimate  this  figure  for  a  typical  forest,  a  value  of  6%  is  assumed  here. 
The  moisture  content  and  salinity  used  axe  representative  of  leaves  at  the  height  of 
the  growth  season.  The  values  calculated  by  the  model  compared  favorably  with 
experimental  data  found  in  the  literature3, 17,20,24,25. 


Frequency  (MHz) 

50 

100 

200 

400 

600 

800 

1300 

2400 

3200 

f' 

44.6 

44.7 

44.6 

44.5 

44.3 

44.3 

44.0 

43.4 

43.0 

€» 

209.0 

105.0 

53.0 

_ 

27.0 

19.0 

15.0 

11.0 

10.0 

11.0 

Table  IB.  ('  and  c"  of  broad  leaves  and  needles  for  leaves  with  moisture  content 
of  65%,  salinity  of  6%,  and  ambient  air  temperature  of  28’C. 
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DENSITY  OF  FOREST  COMPONENTS 


Having  specified  values  for  the  electrical  properties  of  forest  constituents, 
quantitative  values  describing  the  density  of  the  various  foliage  components  must  be 
established.  A  typical  parameter  used  to  define  the  density  of  foliage  is  the  volume 
fraction  that  the  foliage  components  occupy  in  a  typical  forest.  The  volume  fraction  is 
defined  as  the  ratio  of  volume  occupied  by  foliage  components  to  the  total  volume  of 
the  layer  containing  the  forest  environment.  Quantitative  measurements  of  the 
volume  fractions  of  forests  axe  scarce,  and  determination  of  an  estimate  requires  the 
development  of  sophisticated  modeling  programs.  Although  models  are  available  to 
determine  the  volume  fractions  of  forest  stands3,22,24,  actual  values  of  volume  fractions 
can  be  expected  to  vary  greatly  depending  on  the  region,  season  and  weather  pattern. 
The  model  in  this  report  calculates  the  effective  parameters  as  a  function  of  the 
volume  fraction  (Figures  1  through  4)  which  allows  determination  of  £*  based  on  an 
estimate  of  the  volume  fraction  of  a  forest  stand  of  interest.  Therefore,  the  only  task  is 
to  determine  a  suitable  range  of  volume  fractions  over  which  to  do  the  calculation  of  e* 
that  will  be  representative  of  a  majority  of  forest  environments. 

From  the  literature  reviewed,  it  was  found  that  in  the  typical  mature  forest  the 
wood  will  occupy  approximately  1.0%  of  the  total  volume2,3,22,24.  In  some  cases,  well 
managed  forests  can  attain  wood  volume  fractions  in  the  vicinity  of  5.0%.  This  is 
considered  an  extreme  case.  For  the  purpose  of  this  work,  calculations  will  be  done  for 
wood  volume  fractions  over  a  range  of  0.1%  to  5.0%. 

The  volume  fraction  of  broad  leaves  or  needles  was  typically  reported  to  be 
between  10  to  25  percent  of  the  volume  fraction  occupied  by  the  wood.  Therefore,  the 
volume  fractions  of  leaves  in  this  report  will  be  consistently  reported  as  10%  of  the 
total  volume  fraction,  which  is  considered  typical3.  It  is  the  intent  to  use  a  range  of 
volume  fractions  that  is  significantly  broad  so  that  most  forest  stands,  with  similar 
vegetation  composition,  can  be  represented. 
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areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
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